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A dynamical study of the decay of a metastable state by quantum tunneling through an 
anisotropic, non separable, two-dimensional potential barrier is performed by the numerical so- 
lution of the time-dependent Schrodinger equation. Initial quasi-stationary proton states are chosen 
in the framework of a deformed Woods-Saxon single-particle model. The decay of two sets of states 
corresponding to true and quasi level-crossing is studied and the evolution of their decay properties 
as a function of nuclear deformation is calculated around the crossing point. The results show that 
the investigation of the proton decay from metastable states in deformed nuclei can unambiguously 
distinguish between the two types of crossing and determine the structure of the nuclear states 
involved. 



I. INTRODUCTION 



The problem of the decay of a metastable state is of major importance in many different fields like nuclear chemistry, 
diffusion processes, spontaneous radioactivity, etc (for a review, sec UH). A metastable state, or quasi-stationary 
state, is defined as a state of local stability which decays with a rather large but finite lifetime toward a true stable 
minimum. When the temperature of the decaying system is low, quantum tunneling is the dominant decay process. 

The one-dimensional tunneling process has been studied intensively from the very beginning of quantum mechanics 
and it is now quite thoroughly understood, at least within a stationary approach ]|,|||. On the other hand, little is 
known about multidimensional tunneling although some progress has been made through the last two decades [^Q- 
Most of these approaches are restricted to special cases and cannot constitute a full and general framework for the 
study of multidimensional tunneling. 

The recent, huge progress of computer performances makes possible the direct treatment of this complex quantum 
phenomenon through the numerical solution of the Time-Dependent Schrodinger Equation (TDSE) for well prepared 
initial quasi-stationary states, and hence the study of the dynamics of their decay. This procedure has been applied 
successfully to a-decay || and fission (!(]] in a one-dimensional model. A transient period at the very beginning of 
the decay process has been found that corresponds to deviations from the well-known exponential decay law ]Tl[ | . The 
relation between the asymptotic value of the time-dependent decay rate and the commonly used stationary value was 
also pointed out. 

Recently, this method has been generalized to two dimensions to study the tunneling of quasi-stationary states 
through an anisotropic, non separable, potential barrier Jl2[ . Applied to proton emission from excited states in 
deformed nuclei, it has been shown that the escape directions chosen by the proton state depend mainly on its 
quantum numbers, i.e., on its spatial distribution, rather than on the features of the potential. The important role 
played by the distribution of the angular momentum and its variation in time for the determination of the proton 
decay rate have been also pointed out. Along with the experimental observation of nuclei beyond the proton dripline 
jl3| , several stationary calculations for the half-lives of spherical |l4|] and deformed [^5) ground state proton emitters 
have been reported. 

It is well known that as a function of nuclear deformation, some single-particle levels cross each other [jl6). Two 
distinct types of crossing can occur: 1) the quantum numbers (spin projection A and parity 7r) of the crossing levels 
are the same but belong to two different oscillator shells (AN — 2) ; in this case, the levels repel each other and their 
wave functions strongly mix near this quasi-crossing point. 2) the quantum numbers are different ; in this case, the 
two levels do not interact and we deal with a true crossing point [pT[JT^] . 
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Therefore one roughly expects the proton decay properties of the two levels: 1) to be similar for quasi-crossing, the 
only difference coming from the small 'interaction' energy ; 2) to be very different for true crossing, although the 
two levels have exactly the same energy. This implies that the proton decay from excited states in deformed nuclei 
can represent a new tool for studying the level-crossing phenomenon and defining its type. Using the numerical, 
time-dependent approach mentioned above, the present paper offers a detailed study of the decay properties of single 
proton excited states for each of the two cases at and around the crossing deformations. 

In the first part, the numerical formalism used for solving the TDSE in two dimensions is described as well as 
the physical quantities accessible with this approach. The deformed single-particle model used to obtain the proton 
quasi-stationary states inside the nucleus is described. In the second part, the time and deformation dependencies 
of the numerical results are presented and discussed. Finally, a summary, a conclusion and possible developments of 
this work are given. 



II. FORMALISM 



A. Time-dependent approach 

We are interested in describing the evolution of an individual proton interacting with an axially symmetric average 
potential created by the nucleons of the daughter nucleus. The total wavefunction representing the proton can be 
written in cylindrical coordinates as 

V<z,p,<M) = /(z,p,i) e ,;A ^ (i) 

where A is the projection of the total angular momentum on the axis of symmetry (z-axis) and it is a good quantum 
number. Hence, the TDSE becomes 

ih^f{z, P) t) = H(X,p){(t,p,U) (2) 

where the hamiltonian Ti is 



2, 



h 2 ri d d 2 d 2 A 2 



p dp dp 2 dz 2 p 2 
+V pA (z,p). 



(3) 



V p a is the potential representing the interaction between the individual proton and the remaining nucleons of the core 
daughter nucleus. 

In order to follow the time evolution of an initial state ip(f 7 t = 0), the TDSE (||) has been solved numerically within 
a finite spatial grid using a multistep predictor method called MSD2 Jl9| . This method presents a good compromise 
between stability, accuracy and efficiency for this type of numerical problems (linear, stationary potential). Within 
this approach, one has access to the wavefunction ip(f, t) at any time t and then to the following physical quantities : 

• Total tunneling probability given by the fraction of the wavefunction located beyond the top of the two- 
dimensional potential barrier (potential ridge) at time t 

P tun (t) = 2tt / \f(z,p,t)\ 2 pdpdz; (4) 
Jv oui 

• Tunneling angular distribution estimated in spherical coordinates (r, 9, <fi) and defined as 

Ptun(t,e) = ^= \f(r,8,t)\ 2 r 2 dr, (5) 



dn 



r„„ t (9) 



where r out (6) is the radial position of the potential ridge in the direction 8. Note that 



2;r / ^L sin Q dB = Ptun - ( 6 ) 
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• Total decay rate related to the total tunneling probability by 

X(t) = 1 dPtun . ,ys 

1 — Ptun dt 

• Mean value of angular momentum defined as 

< L 2 > (t) = < iP{t)\L 2 \*P{t) > 

~<f>(t)(<I>(t) + l). (8) 

This quantity is important for the determination of the tunneling rate as well as of the angular momentum 
exchange during decay (non-spherical potential). 



B. Quasi-stationary states 

To accomplish the numerical scheme proposed above, knowledge of the initial metastable wavefunction is necessary. 
A quasi-stationary state is defined as an infinitesimally modified eigenstate of the hamiltonian. In a one-dimensional 
model, the two-potential formalism p0| can be used straightforwardly ||. In our two dimensional approach, the 
numerical difference between the preparation of the initial quasi-stationary state as an expansion in an orthonormal 
basis (diagonalization) and the resolution of the TDSE on a discretized spatial grid has been proven to be sufficient 
to create quasi-stationary states without modifying the potential |Q. 

In the case of an axially and reflexion symmetric potential, the quantum numbers labeling each single-particle 
proton state are (Air) with A the projection of the total angular momentum on the axis of symmetry and 7r the parity. 
No spin-orbit term is present in our calculations. 



C. Single- particle potential 

The interaction between the emitted proton and the daughter nucleus is represented by the average potential 
V p a- Obviously, this potential depends on the deformation of the parent nucleus whose shape has been described 
by Cassinian ovals pl| . Only one parameter e is sufficient to describe shapes from spheres (e = 0) to scission point 
configurations (e = 1). 



Hence, V£ A can be written as 



V p 5 A = V £ nucL +V £ CouL . (9) 



The nuclear part is a deformed Woods-Saxon potential 

v E nucL (*, p) = t- — y° , » (io) 

1 + exp(d e (z,p)) 

where d e (z, p) represents an approximation of the distance between the point [z, p) and the deformed nuclear surface, 
Vo the depth and a the diffusivity of the nuclear potential. 

The Coulomb part V e is obtained assuming that the nuclear volume is filled with a uniform charge distribution 
with a sharp edge. 



III. NUMERICAL RESULTS 

As in our earlier work fl^| , the numerical procedure described above has been applied to the hypothetical proton 
emission from excited states in a nucleus characterized by Z — 82 and A — 208. It is worth noting that, qualitatively, 
the results presented are independent of the choice of the decaying nucleus. While in this previous paper we concen- 
trated on the general dependence of proton tunneling on time, nuclear deformation and quantum numbers, here we 
will analyze the behaviour of two selected sets of pair states around the true-crossing and quasi-crossing deformations, 
respectively. 

The parameters used for the potential V p a have been chosen according to Ref . (2^] . The depth Vq has been modified 
in order to obtain states lying below the top of the two-dimensional barrier : Vo — —62.66 MeV for the case of 
true-crossing, and Vo = —58.7 MeV for the case of quasi-crossing. 



3 



A. True-crossing states 

We studied the case of two single-particle states with different quantum numbers, ijj^h (Att = 3—) and ipik(A/K = 
3+)B 

At the deformation etc = 0.25, the two corresponding energy levels cross, i.e., undergo an accidental degeneracy, 
E2h — Eik- This situation arises because of the multidimensionality of the problem. 

In Fig. [l], the square moduli of the two crossing states are plotted in the half-cylindrical plane for three different 
values of e at and around crossing. 

This figure shows clearly that the two wavefunctions don't mix at all during their crossing since they keep their 
spatial distribution constant as the potential deformation is changing. 

Using the procedure described in Sec. ||, we followed the time-evolution of these two initial states for the three different 
deformation values e = 0.23, 0.25, 0.27. To illustrate the numerical results obtained, the square wavefunctions ip2h 
and ipik (at etc) have been plotted in Fig. || at five different times during the calculation, showing the tunneling of 
these quasi-stationary states through the two-dimensional anisotropic barrier. 

In Fig. H are plotted the time evolution of the decay rate as well as the tunneling probability for the two states at 
the three deformations considered. 

Two different stages can be distinguished in the time evolution of A There exists an initial transient time 
corresponding to the "acclimatization" of the quasi-stationary state to its new environment and during which the 
decay rate A undergoes strong oscillations. After that, A is almost constant in time (in the following, this value is 
denoted A(ooV). This second stage corresponds to an exponential decay. 

From Fig. one can extract two important results. First, for each state, the decay rate and the tunneling probability 
evolve very smoothly as the deformation of the potential increases. As expected, the two states don't influence each 
other during crossing. 

Second, at a given deformation e, the two asymptotic decay rates A(oo) differ by about one order of magnitude 
although the two states are degenerate. In fact, this important difference can be explained by the difference in the 
"residual" angular momenta < I > of the two decaying states. Although not plotted here, < I > (t) is almost constant 
in time for both states but at two quite different values : < I >~ 5.26 for ip 2 h and < I >~ 7.56 for ip lk . As shown 
before |Q, this quantity plays a major role in the determination of the decay rate. 

From the knowledge of the wavefunctions in the cylindrical plane, we have inferred the tunneling angular distribu- 
tions Ptun(t, 0) at different times. The results of such calculations are shown in Fig. ^ for the two crossing states at 
four different times and at the deformation etc = 0.25. Obviously, the two angular distributions are really different 
from each other and are strongly correlated to the spatial distributions of the corresponding initial quasi-stationary 
states. Moreover, it is worth noticing that, more generally, an emission at 9 — 90 deg is always related to a state 
symmetric with respect to the p-axis (cf. -02/1 state). Contrarily, if the state is asymmetric (cf. ipik state), the emission 
is hindered at 8 = 90 deg. 

In conclusion, for a true crossing, one can determine which of the two metastable single-proton states is populated. 

B. Quasi-crossing states 

Now we study the influence of the mixing of two orbitals having same A and n but AiV — 2 in the single-proton 
wave function on the tunneling behaviour. 

Without loss of generality, we chose two states with (An) = 3+ whose levels cross at deformation e qc — 0.458. 
Figure [s] represents the square moduli of the corresponding single-particle wave functions at and around the critical 
value e qc , while the evolution of their energy and angular momentum as a function of the deformation parameter e is 
plotted in Fig. ^. At e qc the difference in energy is the smallest (AE qc ~ 290 keV) and the two "residual" angular 
momenta are equal numerically (A < I > qc = O.Of). 

Following the wave function corresponding to the lower (higher) energy E < (£>) vs the deformation e, one notices 
that its spatial distribution evolves very rapidly but continuously. Asymptotically, i.e., for e 3> e qc , they have indeed 
interchanged their spatial distributions. 

The time-dependent decay properties of these states are presented in Figs. |7] and |^ for three different deformations 
e. The shapes of the angular distributions (Fig. ||) for different deformations seem to be again related to the spatial 



1 The indices "2h" and "lk" for the proton states have been chosen using the usual spectroscopic notation according to the 
angular momentum value of the corresponding spherical states (1 = 5 and I = 8, respectively). 
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distribution of the corresponding quasi-stationary states. At the point e qc , the tunneling angular distributions for the 
states tend to coincide. 

The variation of the decay rate from one deformation to another (Fig. |^) is clearly not correlated with the corre- 
sponding variation of the energy. However, the ratio A^ > (oo)/Xe < (oo) (of about 3 at e qc ) can be explained by the 
difference in energy {E > — E < = 290 keV) using stationary one-dimensional WKB estimates. 

To study in more detail the dependence of A on nuclear deformation, time-dependent calculations have been per- 
formed for the same quasi-stationary states at several deformations on both sides of e qc . The asymptotic values of 
the decay rates A(oo) have been plotted in Fig. || as a function of the deformation e. It is clearly shown that at the 
quasi-crossing point (e = e qc ), the two asymptotic decay rates are well separated. Contrarily, at e = 0.4, the two 
A(oo) are almost equal. This effect can be understood in terms of the competition between the energy of the decaying 
state and its "residual" angular momentum value < I >. For e < e qc , the higher energy state has the higher < I > 
value. On the other hand, for e > e qc , the higher energy state corresponds to the smaller < I > value. Both effects 
act in the same direction, strongly enhancing the decay rate. 

In conclusion, for a quasi-crossing, one can determine the energy splitting and whether the proton has jumped into 
the upper level or not (Landau-Zcner effect). 



IV. SUMMARY 



The two-dimensional Time-Dependent Schrodinger Equation has been solved numerically for prepared initial single- 
particle quasi-stationary states decaying through an anisotropic, non separable, potential barrier. The special cases of 
true crossing and quasi-crossing have been studied within this approach. Time-dependence of quantities like tunneling 
probability, decay rate and tunneling angular distributions have been calculated. 

It has been shown that two (accidentally) degenerate proton states decay very differently due to their different spatial 
distribution and "residual" angular momentum. Hence, one-dimensional semi-classical approximations of their decay 
rate are not expected to be reliable. 

The mixing of two orbitals having the same main quantum numbers in the single-proton wave function strongly 
influences the decay behaviour. At the quasi-crossing deformation, the difference between the proton decay rates from 
the two states involved can be entirely explained by the difference in their energy. A nearby deformation, where these 
two rates are the same, can always be found due to the combined effects of the residual angular momentum and of 
the energy of the decaying states. 

Proton decay in strongly deformed nuclei has been recently observed [p3| . Extensions of this type of work would 
allow a detailed study ot the level-crossing phenomenon, including the structure of the nuclear states involved. 

Further developments of the present work should include generalization to proton decay through a time-dependent 
potential barrier since only a dynamical approach could tackle such a problem. There are indeed cases in which the 
parent and daughter nuclei have very different deformations, e.g., prolate-oblate transitions or emission during nuclear 
fission. 
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FIG. 1. Square modulus of the two crossing states at and around the crossing deformation etc = 0.25. Dashed lines represent 
the deformed ridge of the potential for each deformation e. 
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FIG. 2. Time evolution of \1p2hf and of \ipikf at the deformation e tc = 0.25. The time t is given in units of 10 -22 sec. For 
clarity, contour lines with values greater than 4 x 1CF 5 (3 x 1CF 6 ) in the upper (lower) row are not plotted. The potential ridge 
is drawn with dashed lines. 
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FIG. 3. Time-dependent decay rates A and tunneling probabilities Ptun for the two crossing states ip2h and ipik obtained for 
three nuclear deformations. The asymptotic decay rates A(oo) correspond to etc = 0.25. 
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FIG. 4. Time-dependent tunneling angular distributions for the two crossing states ip2h and ipik- The time t is in units of 
1(T 22 sec. 
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FIG. 5. Contour plots of the two mixing wavefunctions at and around e qc = 0.458. Dashed lines are the deformed ridges. 
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FIG. 6. Average energy and "residual" angular momentum of the two quasi-crossing states as a function of deformation. 
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FIG. 7. Time-dependent tunneling probabilities Ptun and decay rates A for the two quasi-crossing states at three deformations. 
The values of the corresponding "residual" angular momenta are indicated in the lower part. 
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FIG. 8. Time evolution of the tunneling angular distributions dPtun/dQ. for the two quasi-crossing states at three deforma- 
tions. The distributions have been computed at four times : t = 5, 10, 15 and 20 x 10 -22 sec. 
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FIG. 9. Asymptotic decay rates A(co) for the two initial (3+) states vs the deformation e. 
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